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Abstract: 


of a point given in geographic. military, state planar. or Universal Transverse 
Mercator form to coordinates represented in a second form chosen from these 
four coordinates types. 


This set of twelve FORTRAN subroutines will convert the coordinates 
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1.0 INTRODUCTION 


In the Installation Restoration (IR) program measurements of variables such as depth 
to top of interval in a coring must be related to a specific location on the earth’s surface. This 
is accomplished through the use of coordinates: geographic, state planar, Universal Transverse 
Mercator, and military. Each installation will use that coordinate system which is most con- 
venient; therefore, the programmer will be concerned with conversion from one coordinate svs- 
tem to another, This manual will show you how to use the subroutine package. If you under- 
stand the above coordinate systems. you may skip section 2. 


Note: Your comments on this manual can help improve it: please direct all comments to 
Potomac Research, Inc., at our Edgewood Area office. 
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2. LOCATING POINTS 


Ww 


0 General 


The earth’s surface is considered to result from the rotation of an ellipse about the 
north-south polar axis. This solid is called the spheroid. The Clarke 1866 spheroid is cur- 
rently used for the continental United States. Dimensions are shown below: 


NP 
b a = semimajor axis 
a 6.378.206.4 meters 
b = semiminor axis 
6,356,583.8 meters 
SP 


Figure 1. Clarke 1866 Spheroid 


A point such as position A in Figure 2 is located by describing it in terms of two numbers called 
coordinates. 


Geographic coordinates can be related to the spheroid alone by means of the formulae of 
geodesy. For the other coordinate systems, the point A has been transferred from the spheroid 
toa map projection. The map projection is a flattening of the surface upon which the spheroid 
is projected according to the projection type. We desire a flat surface because it is easier to work 
with than a curved surface. - 


To locate that point A in Figure 2, we measure its signed distances from two lines called 
axes and record the numbers in order as the coordinates of A. Our x and y axes will always inter- 
sect at an angle of 909, 


X axis 


+x 


coordinates + X, + Y 


Y Figure 2. General Coordinate System 
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2.1 Geographic Coordinates 
ae 
On the spheroid the axes are curved lines. The x axis is the circle equidistant from 
the North and South Poles named the Equator. The v axis is that half-ellipse through the 
North and South Poles that passes through Greenwich. England. named the Prime Meridian, 
Distances are meusured in degrees, minutes, and seconds of arc. We shall measure from 0 
to 90° north or south of the Equator and shall call this coordinate latitude. We shall mea- 
sure from 0 to 1SO0° west or east of the Prime Meridian and shall call this coordinate longi- 
tude. The precision of our measurements shall be to the nearest 0.01 second. 
Geographic coordinates cf point A are written like the following example: 
399 23° 51.37"N 76° 18’ 06.63” W 
Sometimes. you will see these symbols: 
@, latitude 


A, longitude 


Also, geographic coordinates may be called geodetic. 


2.2 State Planar Coordinates 
Two types of projections are used in the state planar coordinate system: 


Lambert projection — ‘‘wide” state such as Connecticut 


Transverse Mercator projection — “tall”? state such as Vermont 
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nay Ire e the ne projection to caver a state. as the width is limited to 1F 
it may require more than one projection to caver a state. as the widti lintited to TEs 


coded us four-digit numbers in subroutine FLALPA and ZONES. The y ann fs the central 
meridian of the projection: it is assigned a large positive number so that no coordivate will 
have a negative value, Its equation is x = p. The x anis is a line perpendicular to thes 

aXis and it is sufficiently south of the zone se that no y coordinate will have a newative value. 
Its equation is y = O. We shall meusure distance in feet from the ¥ axis te point Az cull 


this distance X'. X coordinate is 


X =p + X! 


Note: Use + if point A is east of the central meridian: use - if point is west. 


We shall measure distance in feet from the x axis to point A and shall call this coordinate 


Y. 
State planar coordinates of point A are written like the following example: 
Zone 4126 X = 997,307 Y = 570,443 
X=p 
vy =O 

central meridian 

Figure 3. State Planar Coordinate System 
| 
2.3 Universal Transverse Mercator Coordinates 


The Universal Transverse Mercator (UTM) projection is used from 80°S through 
84°N latitude. The earth is divided into 60 zones 6° wide starting at 180°W longitude and ' 
proceeding east. Zones are numbered from | to 60 as in Figure 4. 
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180°W 174° 168° eo. 78OW 729 . 168° 1749 180P°E 


Figure 4+. UTM Zones 


The y axis of each zone is its central meridian: it is assigned a value of $00,000 meters. Its 
equation is E = 500,000. 


The x axis of each zone is the Equator; in the northem hemisphere, it is assigned a value of O 
meters; in the southern hemisphere, 10,000,000 meters. Its equation is N = 0 (northern); 
= 10,000,000 (southern), We shall measure distance in meters from the y axis to point A: 
call this distance E', The false easting coordinate is 
E = 500,000 + E' 


Note: Use + if point A is east of central meridian; use - if point A is west. 


We shall measure distance in meters from the x axis to point A and shall call this northing 
distance N', The false northing coordinate is 


N =N! (northern hemisphere) 
or N = 10,000,000 - N' (southern hemisphere) 


UTM coordinates of point A are written like the following example: 
Zone +18 E = 387,900 N =4,361,500 


Note: + zone means northern hemisphere 
- zone means southern hemisphere 
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Military coordinates are based on the OPM svstem described in section 2.3, 


mone is divided into S© high areas called) Guodrangles except from TION to S4ON where 


the quadrangle is 6 ~ TDS The row dimension of a quadrangle is lettered from C through 


XN comitting Pand Oi as in Figure 3.0 A quadrangle is identined by zone number und row, 


S4ON 
720 x 
64° W 
569 x 
4§0 U 
40° T 
320 S 
240 R 
16° Q 
go P 
0° N 


64° . 
720 D 
80°S Cc 


Figure §, Quadrangles of a UTM Zone 


These quadrangles are then divided into squares which have 100,000 meters on each side. 

A 100.000-meter square is identified by two letters: x designator, y designator, Compute 

r = MOD (zone, 3)*. If r equals 0, then the x designator is lettered S through Z proceeding 
“ 


east: if r equals 1, then it is lettered A through Hy if r equals 2, then it is lettered J through 
R (omitting O). 


If the UTM zone is odd, then the y designator is lettered L through V, A through K 
(omitting I and O) starting at a multiple of 2,000,000 meters and proceeding north; if the 


UTM zone is even, then it is lettered R through V, A through Q (omitting I and O) as in 
Figure 6. 


*Fortran MOD function. 


POTOMAC RESEARCH, INCORPORATED 


ies 


SOO 001) Meters 


Row § 
(from 320N to 40UN) 


—— 4.000.G00 Meters 


329N 


750W 


Zone 18 
Note: Because of convergence of the meridians. letters S and Z have been cut off. 


Figure 6. 100,000-Meter Squares of an Even Zone 


Since the point A hus been located within 100.080 meters by means of zone number, quadrangle row, 
and 100.000-meter square, We use only part of the casting and nortiing coordinates: viz.. the mehtmost 
tive digits of each. 


Military coordinates of point A are woitten as a character string using the elements discussed above 
trom lett to right like the following example: 


18SUU8790061500 
For readability, the coordinates above may be sepurated: 
18 SUU 8790061500 


Note: + UTM zone means northern hemisphere 
- UTM zone means southern hemisphere 


This sign is returned as a bonus from the appropriate subroutine: it is not needed except in level 2 
and 3 subroutines. 
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3. USER PROCEDURES 


3.0 General 


For each type of coordinate conversion, one subroutine has been provided so that the subroutine 
call will be shorter than that arising from one general purpose subroutine, 


a Use this UNIVAC control statement: 
LIB) = =IR*#CVCOORD. (with &© MAP) 
b. Name the subroutine which you desire to call: 
COORxy 
where COORxy will convert from type x to y 
x, y is one of the letters G,S, M, U 
(G = geographic, S = state planar, M = military, U = UTM). 
c. Check the calling sequence in sections 3.1 through 3.12. Be sure that your input data con- 
form to the formats shown therein: e.g.*, longitude LAM is one double precision number formed from 


degrees, minutes, and seconds: 


76°18°6.63”"W 0761806.63D0 


(longitude format = dddmmss.ssDO) 
d. Use Table I to get your installation input code (INST). 


e. After each subroutine call, test the error code (JERR) returned. If JERR equals 0, then the 
conversion hus been successful: if JERR does not equal 0, then use Table III to find your error. 


* Aliso, these variables must be in the correct format: latitude PHI (similar to LAM), quadrangle row and 
100.000-meter square QSQ (3 letters), military coordinate string XYMIL (10 digits D.P.). 
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Table f. Installation Codes 


Code Installation, 

| AL Alabama Army Ammunition Plant 

| BA Badger Army Ammunition Plant 
EA : Edgewood Area 
FA Frankford Arsenal 
HA Hawthorne Naval Ammunition Depot 
LS Lone Star Army Ammunition Plant 
PB Pine Bluff Arsenal 

| RA Redstone Arsenal 
RM Rocky Mountain Arsenal 
vO : Volunteer Army Ammunitions Plant 
ws Weldon Spring Chemical Plant 


Table II. Zone Numbers{provided in subroutine ZONE2) 


Installation State Zone . UIM Zone 
AL 3101 +16 
BA 5751 +16 
EA 4126 +18 
FA 5151 +18 
HA 4651 +11 
LS 5351 +15 
PB 3251 +15 
RA 3101 +16 
RM : 3451 +13 
vo 5301 +16 
WS 4401 +15 
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Note. T= input variable. QO = output variable 


3.1 Geographic to Military Coordinates. 


CALL COORGM (NST. PHI, LAM. UZONE, OSQ, NYMIE, JERR) 


INST -- YJTHOLLERITH INSTALLATION (Table 1. page 9) 
PHI -- ID.P. LATITUDE 
LAM -- ID.P. LONGITUDE 
UZONE -- OINTEGER +OR — UTM ZONE NUMBER 
QSO -- OHOLLERITH QUADRANGLE ROW & 100.000-M. SQUARE 
AYMIL -- OD.P. COORDINATE STRING, MIL. 
JERR -- QINTEGER ERROR CODE (Table HI. page 15) 
Sr Geogruphic to State Planar Coordinates 


CALL COORGS (INST. PH], LAM, SZONE, XSTP, YSTP, JERR) 


INST -- JHOLLERITH INSTALLATION (Table I, page 9) 
PHI -- ID.P. LATITUDE 

LAM -- ID.P. LONGITUDE 

SZONE -- QINTEGER STATE ZONE NUMBER 

XSTP + -- OD.P. X COORDINATE, STATE 

YSTP -- OD.P. Y COORDINATE. STATE 

JERR -- OINTEGER ERROR CODE (Table II, page 15) 


3.3 Geographic to UTM Coordinates, 


CALL COORGU (INST, PHI, LAM, UZONE, EAST, NORTH, JERR) 


INST -- THOLLERITH INSTALLATION (Table I, page 9) 
PHI -- ID.P. LATITUDE 

LAM -- ID.P. LONGITUDE 

UZONE -- OINTEGER +OR —UTM ZONE NUMBER 
EAST -- OD.P. EASTING 

NORTH -- OD.P: NORTHING 

JERR -- OINTEGER ERROR CODE (Table II], page 15) 
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3.4 Military to Geographic Coordinates, 


CALL COORMG UNST, OSQ. XYMIEL. HAT PHL H2. LAM. IP RR) 


INST 
OSO 
NYMIL 


HH] 
PHI 
Hi 
LAM 
JERR 


ne 
wry 


THOLLERITH INSTALLATION (Tobie To pare 9) 
THOLLFRITH QUADRANGLE ROW & 10U.000-M. SOU ARE 
ID. P. COORDINATE STRING, MIL. 


OHOLLERITH HEMISPHERE (N or S) 


OD.P. LATITUDE 
OHOLLERITH HEMISPHERE (Wor F) 
ODP. LONGITUDE 


O INTEGER ERROR CODE (Table HL. page 15) 


Military to State Planar Coordinates. 


CALL COORMS UNST. QSQ, XYMIL. SZONE, XSTP, YSTP. JERR) 


INST 
QSQ 
XYMIL 


SZONE 
XSTP 
YSTP 
JERR 


IHOLLERITH INSTALLATION (Table I. page 9) 
IHOLLERITH QUADRANGLE ROW & 100.000-M. SQUARE 
ID. P. COORDINATE STRING. MIL. 


O INTEGER STATE ZONE NUMBER 

0D.P. X COORDINATE, STATE 

OD.P. Y COORDINATE, STATE 

0 INTEGER ERROR CODE (Table II, page 15) 


3.6 Military to UTM Coordinates. 


CALL COORMU (INST, QSQ, XYMIL, UZONE, EAST, NORTH, JERR) 


INST 
QSQ 
XYMIL 


THOLLERITH INSTALLATION (Table I. page 9) 
IHOLLERITH QUADRANGLE ROW & 100,000-M. SQUARE 
ID.P. COORDINATE STRING, MIL. 


0 INTEGER + OR — UTM ZONE NUMBER 
0D. P. EASTING 

0D.P. NORTHING 

0 INTEGER ERROR CODE (Table III, page 15) 
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3.7 State Planar to Geographic Coordinates. 
CALL COORSG (INST. NSTP. YSTP. HI. PHI HE. LAM. JERR) 
INST - THOLLERITH INSTALLATION (Table 1. page 9) 
NSTP - 1D. P. X COORDINATE, STATE 
YSTP - ID.P. Y COORDINATE, STATE 
Hl - OHOLLERITH HEMISPHERE (N or S) 
PHI - OD.P. LATITUDE 
H2 ; OHOLLERITH HEMISPHERE (W or E) 
LAM - OD.P. LONGITUDE 
JERR - O INTEGER ERROR CODE (Table III. page 15) 
38 State Planar to Military Coordinates. 


CALL COORSM (INST, XSTP, YSTP, UZONE, QSQ, XYMIL. JERR) 


INST - 
XSTP - 
YSTP - 


UZONE - 
QSQ 2 
XYMIL - 
JERR - 


3.9 State Planar to UTM Coordinates. 


] HOLLERITH 
1D. P. 
1D. P. 


O INTEGER 

0 HOLLERITH 
OD.P. 

O INTEGER 


INSTALLATION (Table I. page 9) 
X COORDINATE, STATE 
Y COORDINATE, STATE 


+ OR —UTM ZONE NUMBER 
QUADRANGLE ROW & 100,000-M. SQUARE 
COORDINATE STRING, MIL. 

ERROR CODE (Table HII, page 15) 


CALL COORSU (INST. XSTP, YSTP, UZONE, EAST, NORTH, JERR) 


INST - 
XSTP - 
YSTP - 


: yA. = a I 


1 HOLLERITH 
iD. P. 
ID. P. 


0 INTEGER 
0D.P. 
0D. P. 
O INTEGER 


INSTALLATION (Table I, page 9) 
X COORDINATE, STATE 
Y COORDINATE, STATE 


+ OR — UTM ZONE NUMBER 
EASTING 

NORTHING 

ERROR CODE (Table III, page 15 ) 
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3.10 UTM to Geographic Coordinates. 


CALL COORUG UNST. EAST. NORTH. HL. PHE HO. DAM. JE RRG 


INST -- THOLLERITH INSTALLATION (Table [. pege 9) 
EAST -- ID.P. EASTING : 
NORTH -- ID.P. NORTHING 

mW -- QOHOLLERITH HEMISPHERE (N or S) 

PHI = ODP: LATITUDE 

H2 -- QHOLLERITH HEMISPHERE (W or F) 

LAM -- OD.P. LONGITUDE 

JERR -- OINTEGER ERROR CODE (Table HI. page 15) 


3.11 UTM to Military Coordinates. 


CALL COORUM (INST, EAST. NORTH, UZONE, OSQ. XYMIL,. JERR) 


INST -- YTHOLLERITH INSTALLATION (Table I, page 9) 

EAST -- ID.P. EASTING 

NORTH -- —D.P. NORTHING 
UZONE -- OINTEGER +OR-UTM ZONE NUMBER 

QSQ -- QOHOLLERITH QUADRANGLE ROW & 100,000-M. SQUARE 
XYMIL -- OD.P. COORDINATE STRING, MIL. 

JERR --  OINTEGER ERROR CODE (Table HI. page 15) 


3.12 UTM to State Planar Coordinates 


CALL COORUS (INST, EAST, NORTH, SZONE. XSTP, YSTP, JERR) 


INST -- TPHOLLERITH INSTALLATION (Table I, page 9) 
EAST -- ID.P. EASTING 

NORTH -- ID.P. NORTHING 

SZONE -- OINTEGER STATE ZONE NUMBER 

XSTP -- OD.P. X COORDINATE, STATE 

YSTP -- OD.P. Y COORDINATE, STATE 

JERR -- OINTEGER ERROR CODE (Table III, page 15) 
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3.13 Example. 


You have an input file of well focations in state plenur coordinates from Rocky Mountain Arsenal, 
The following program fragment will convert them te UTM coordinates: | 
INTEGER UZONE 
DOUBLE PRECISION E.N, XSTP. YSTP | 


200 READ (5.3177, END=9000) XSTP, YSTP 
WRITE (6.3105) XSTP, YSTP 
CALL COORSU (2HRM. XNSTP, YSTP. UZONE. EL N. JERR) 
IF (JERR. EQ. 0) WRITE (6.3105) ELN 
GOTO 200 


9000 CONTINUE 
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Table HI. Error Codes Returned by Subroutine 
Error Code In 
_UERRY Explanation Subroutine: 
0 No errors! COOR__ 
| INVALID INSTALLATION COOR__ 
3 MISSING INPUT COORDINATES COOR __ 
3 MISSING 3 LETTERS -- QUADRANGLE ROW COOR__ 
& 100.000-M. SQUARE 
‘ 
10 UTM ZONE NUMBER > 22 WITH CLARKE 1866 CONVMU 
SPHEROID 
11 INVALID ROW LETTER OF 6° X 89° QUAD- CONVMU 
RANGLE (MUST BE N THROUGH X) 
}2 INVALID X DESIGNATOR OF 100,000-M. CONVMU 
SQUARE (SEE PAGE 6) 
13 INVALID Y DESIGNATOR OF 100,000-M. CONVMU 
SQUARE (SEE PAGE 6) 
14 INVALID NORTHING WITH CLARKE 1866 CONVMU 
SPHEROID 
15 LATITUDE > 849 TMFWD 
16 LONGITUDE > 0.16 RADIANS TMFWD 
17 GRID DISTANCE FROM CENTRAL MERIDIAN TMINV 


> 0.2a WHERE a = semimajor axis of spheroid 
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Table HE. Frror Codes Returned by Subroutine — Continued 


18 RECTIFYING LATITUDE > 1.47 RADIANS TMINY 
19 INVALID STATE ZONE NUMBER CONVGS 
20 LATITUDE > 80° 30° UCPMPEW 
il NORTHING > 8,942,900 M. UTMIN'S 
| ' 
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The twelve subroutines were tested separately on one point within each installation, 


Results were then compared with the coerdinates measured from the map for that installation. 


Results were in agreement. 


Additionally. eycles of forward and backward conversions such as G to §.S to G. ete.. were tested. 
Results agreed within accuracies of 0.0005 meter. 1.2 feet. or 0.02 seconds. 


TEST POINTS 


(Coordinates are listed in the order geographic. state. UTM, military.) 


Aberdeen Proving Ground (Edgewood Area) 


‘ 399 23° 51.37°N 
State zone 4126 
UTM Zone + 18 
+18 SUU 


76° 18° 6.63°W 


X =997,307 
E= 387,900 
8790061500 


Badger Army Ammunition Plant 


439 22° 0.08"N 
State zone 575] 
UTM Zone +16 


+16 TBD 


Franktord Arsenal 
40° 1° 3.317°N 
State zone 5151 
UTM Zone +18 

+18 TVV 


Pine Blutf Arsenal 
34° 19’ 0"N 
State zone 3251 
UTM Zone +15 

+15 SWH 


899 45° 7.92"W 


X= 2,065.894 
E = 277,000 
7700005000 


75° 2° 48,747°W 


X = 2,757,015 
E = 496,000 
9600029500 


92067 0"W 
X = 1,969,804 
E = 582,806 
8280697439 


Y = 570,443 
N= 4,361,500 

Y = 498,207 

N= 4,805,000 

Y = 260,846 


N = 4,429,500 


Y = 600,410 
N = 3,797,439 
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Redstone Anenal 
34° 40° O"N 86° 38° OTW 
State zone 310) X= 259.433 
UTM Zone +10 E 3 
+16 SEP 3359535943 

Rocky Mountain Arsenal 
399 47° 54,99"N 1049 54° 10.98" W 
State zone 3451 X = 2,167,734 Y = 170,047 
UTM Zone +13 E = 508,300 N = 4,405,200 
+13 SEQ 0830005200 

Weldon Spring Chemical Plant 
38° 41° 49.40"N 90° 43° 44,49°—=W 
State zone 4401 X= 434,636 Y = 1,042,715 
UTM Zone +15 E= 697,500 N = 4,285,400 
+15SXN 9750085400 

Alabama Army Ammunition Plant 
33° 19° 59.99” N 86° 1970” W 
State zone 3101 X= 352,369 _ Y= 1,031,023 
UT zone +16 E= 563,595 N = 3,688,257 
+16 SEM 6359588257 

Hawthorne Naval Ammunition Depot 
38° 34° 9.13" N 118° 36’ 49.97°W 
State zone 4651 X = 491,266 Y = 1,390,319 
UTM zone +11 E= 359,400 N = 4,270,000 
+11 SLN 5940070000 

Lone Star Army Ammunition Plant 
33° 26’ 0” N 94° 1370” W 
State zone 5351 X = 3,001,488 Y = 658,370 
UTM zone +15 E = 386.897 N = 3,699,797 
+15 SUG 8689799797 
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Volunteer Arms Ammunition Plant 


BSP 5° 10,907 N SSO 8 34.417 W 


State zone SA01 X= 2.256.446 Y = 253.889 
UTM zone +16) E = 669.300 N = 3.884.000 
+16 SFP 6930084000 
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bs the U.S. Geolegtead Siirves. 


PROGRAMMING NOTES 


General 


“hy roby . roy tp ty tthe wy tue ae ns .oune r Nesp gd od N : a tas 
The twelve basic subroutines use subroutines ddapted from program Dods4 acd JS. tearitten 


Subroutine CONVME was developed bs Potomie Rossurci. Po corpenited. 


Level 3: 


TMCOF 
TMINV 
TMFWD 


Schroutines are organived tato three hicrarchieal levels as follows: 
Level |: 
COORGM COORGS COORGU COORMG 
COORMS COORMU COORSG COORSM 
COORSU COORUG COORUM COORUS 
Level 2: 
CONVGL CONVGS CONVMU ZONI-2 


FLALFA STOD OBMER 
FILA UIMFW OBMINYV 
FILLA UTMIVS GTPLAM 
FILLAA LAMFRD PTGLAM 
FILLBB LAMINV 


Subroutines with the name COOR -- perform the conversion in one direction: subroutines with the name 


CONV-- in two directions. 


Subroutine functions are as follows: 


Subroutine Name 


Function 


COORGM Convert geographic to military coord. 
COORGS Convert geographic to state coord. 

COORGU Convert geographic to UTM coord. 

COORMG Convert military to geographic coord. 
COORMS Convert military to state coord. 

COORMU Convert military to UTM coord. 

POTOMAC RESEARCH, INCORPORATED 
a re ene eine: Sees seh Sees 


COMING, 


COCERSSS 


COORSU 


COOR LG 


COORUM 


COORUS 


CONVGS 


CONVGU 


CONV MU 


GTPLAM 


PTGLAM 


| FLALPA 


FILA 


FILLA 


FILLAA 


FILLBB 


LAMFERD 


LAMINY 


OBMER 


OBMINYV 


STOD 


TMCOF 


TMFWD 


TMINV 


ere ae, 
bed 
(. i ‘ 
€ xt Cee aes . 
Co Mstate RP UTM coord 
Conmecre LEM wesw reer hig ote re 


Consert UTM to military coord. 

Convert UTM to stote ceord. 

Convert geegrephic to state coord. and vice versa 
Convert geourephic te CPM coord. and vice versa 
Convert military to UTM coord. and vice versa 

Convert geographic to state coord. (Michigan Lambert) 
Convert state to geographic coord. (Michigan Lambert) 
Ful ALPHA (array with state cone (nto 


Fill A() ) array — 2 dim 


to 


Fill A () ) array — 2 dim 


Fill AC) ) array — 3 dim 

Fill B¢ ) array — 3 dim 

Convert geographic to state coord. (Lambert) 

Convert state to geographic coord. (Lambert) 

Convert geographic to state coord. (Oblique Mereater) 
Convert state to geographic coord. (Oblique Mercator) 
Convert seconds to deg-min-sec 

Compute coefficients in A(€  ) array 


Convert geographic to UTM coord. 


Convert UTM to geographic coord. 
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Piste 
Subroutine Name. Function 
UCIMEW Convert geograpide to state coord. (trans. Mercator) 
UTMIVS Convert state to geographic coord. (trans, Mercator) 
ZONE? Find state and UTM zones. hemisphere indicators. 


Conventions adhered to within these subroutines: 
Clarke 1866 spheroid 
+ LATITUDE NORTHERN HEMISPHERE 
- LATITUDE SOUTHERN HEMISPHERE 


+ LONGITUDE WESTERN HEMISPHERE 
LOSGITUDE EASTERN HEMISPHERE 


UTM ZONE NOPTHERN HEMISPHERE 
- UTM ZONE SOUTHERN HEMISPHERE 


+ 


5.1 Modifying the Subroutines. 
(* means a moderate programming effort is required.) 
a. To add another installation: Add installation code to IALF (_ ), hemisphere indicators 


to HSIGNI ( ) and HSIGN2 (_), state zone number from subroutine FLALFA to ISZ (_), and the 
UTM zone to 1UZ (+). The above arrays are contained in subroutine ZONE?. 


b. To add aspheroid*: Change each subroutine COORxy so that it will have a new argument, 


JSPHER (spheroid type). Use the spheroid types contained in subroutine CONVGU. Subroutine 
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CONVGS will have to be chanved to recognize spherotd type. Apparently. numerical information 
coneerning the spheroid is stored in arrays called A (2 dim). A ¢3 dim). and Bs dim). 


Array In Labeled Common In Subroutine, 
A (2 dim) CONST FILA. FILLA 
A(3 dim) OMF FILLAA 
B (3 dim) OMI FILLBB 
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